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We consider the structure of the magnetic fields inside the neutron stars in Eddington-inspired Bom-lnfeld 
(EiBl) gravity. In order to construct the magnetic fields, we derive the relativistic Grad-Shafranov equation in 
EiBl, and numerically determine the magnetic distribution in such a way that the interior magnetic fields should 
be connected to the exterior distribution. Then, we find that the magnetic distribution inside the neutron stars 
in EiBl is qualitatively similar to that in general relativity, where the deviation of magnetic distribution in EiBl 
from that in general relativity is almost comparable to uncertainty due to the equation of state (EOS) for the 
neutron star matter. However, we also find that the magnetic fields in the crust region are almost independent 
of the coupling constant in EiBl, which suggests a possibility to obtain the information about the crust EOS 
independently of the gravitational theory via the observations of the phenomena associated with the crust region. 

In any case, since the imprint of EiBl gravity on the magnetic fields is weak, the magnetic fields could be a poor 
probe of gravitational theories, considering many magnetic uncertainties. 

PACS numbers: 04.40.Dg, 04.50.Kd, 04.40.Nr 


I. INTRODUCTION 


The magnetic field is one of the principal properties in the phenomena of the astrophysical objects. In fact, it is believed that 
the magnetic fields can play an important role during supernova explosions, gamma-ray bursts, jets from active galactic nuclei, 
and so on. The existence of strongly magnetized neutron stars, the so-called magnetars, is also suggested via the measurements 
of spin period and its down rate of the central objects in soft gamma repeaters (SGRs) and anomalous X-ray pulsars. According 
to the magnetic dipole model, the strength of surface magnetic fields of the magnetars is considered to be as large as — 10^® 
G QH. From the SGRs, sporadic radiations of 7- and X-rays are observed, while fierce flare activities called the giant flares 
are also detected on rare occasions. In particular, the quasi-periodic oscillations (QPOs) discovered in the afterglow of the giant 
flares give us the evidences of the oscillations of the magnetized neutron stars ||3l. To theoretically explain the QPO frequencies, 
there are many attempts in term of the crustal oscillations 10-0 and/or magnetic oscillations ll9l- [l0 . In any case, in addition 
to the EOS for neutron star matter, the structure of magnetic fields inside the neutron star must be crucial to understand such 
phenomena. 

On the other hand, the gravitational theory must be also imperative to discuss the relativistic objects. The general relativity 
is mathematically beautiful theory of gravity, and its validity has been probed by a lot of experiments and astronomical ob¬ 
servations. However, most of the verifications of general relativity have been done in a weak field regime, such as the Solar 
System 10, while the tests in a strong field regime are very poor. Perhaps, the gravitational theory describing the astronom¬ 
ical phenomena in a strong field regime might be different from general relativity. This is a reason why modified theories of 
gravity are proposed. Since the observable properties could depend on the gravitational theory, one would see the imprint of the 
gravitational theory as an inverse problem 1110. In fact, the science technology is developing increasingly, which will enable 
us to observe the relativistic objects and phenomena around such objects with high precision. Probably, the gravitational waves 
radiated from such a system are also one of them. Through these observations, it is possible to probe the gravitational theory 
00. 

As a modified theory of gravity, EiBl is recently drawing attention in the context of the avoidance of the big bang singularity 
This theory is originally proposed by Banados and Eerreira based on the gravitational action proposed by Ed¬ 
dington ll^ and on the nonlinear electrodynamics by Born and Infeld ll23l] . EiBl is developed according to a Paratini approach, 
where the connection is considered as an independent field, because the field equations contain ghosts in the metric approach 
(|20. The deviation of EiBl from general relativity can be seen only when the matter exists, i.e., EiBl in vacuum is completely 
equivalent to general relativity, and the deviation becomes significant in high density region. Thus, the compact objects are good 
candidates to see such a deviation. Up to now, the compact objects in EiBl are discussed on several occasions and shown the de¬ 
viation in stellar properties from the expectations in general relativity 1201^ . Maybe, via the direct observations of such stellar 
properties, one would distinguish EiBl from general relativity. We remark that in EiBl the curvature singularity can appear at 
the stellar surface for polytropic EOSs ll^ . which must be a problem to solve even thought this theory is attractive. 
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However, in spite of the importance of magnetic effects in astronomical phenomena, the magnetic fields on the neutron stars in 
EiBI have not been considered. There are solely the considerations on the electrically charged black holes in EiBI 
Thus, in this paper, we consider the magnetized relativistic stellar models in EiBI. This discussion could become a fist step 
to examine the phenomena associated with the neutron stars in EiBI. Actually, there are many uncertainties in the magnetic 
properties, such as its geometry and the currents supporting it, even for a given fixed EOS. So, it must be quite difficult to see the 
imprints of the gravitational theory in the magnetic properties, if neutron stars would have different magnetic geometry and/or 
crust properties irrespectively of the theory of gravity. In this paper, to see how the magnetic fields depend on the gravitational 
theory, we especially focus on the axisymmetric dipole configuration of magnetic fields, because such configuration could be 
dominant in the old neutron stars. The additional factors to determine the magnetic properties should be taken into account, but 
we neglect such effects here to simplify the problem. In this paper, we adopt geometric units, c = G = 1, where c and G denote 
the speed of light and the gravitational constant, respectively, and the metric signature is (—, +, +, +). 


II. MAGNETIZED STELLAR MODELS IN EIBI 


Before considering the stellar models in EiBI, we briery mention EiBI. This gravitational theory is obtained from the action 
S given by 


s = y - A\/^) + S'M[g, ^'m], (1) 

where g and denote the determinants of the physical metric 5 ^ 1 , and is the Ricci tensor 

constructed from the connection T^^^, and denotes the matter action depending on the metric g^^, and matter field T'm- That 
is, the matter field is assumed to minimally couple to the metric tensor, i.e., the matter action depends on g^^i, independently 
of the connection T. This theory also has a dimensionless constant A and the Eddington parameter k, which are related to the 
cosmological constant as A = (A — l)//c. Since we especially focus on the asymptotic flat solutions (A = 0) in this paper, 
hereafter we take A = 1. On the other hand, k is constrained from the observations in solar system, big bang nucleosynthesis, 
and the existence of neutron stars l2.5l 1 3^ [361] . The existence of neutron stars can give us the strong constraint on k, i.e., 
|«;| < 1 m® kg“^ s“^ im. Recently, the possibility to constrain k with the terrestrial measurements of the neutron skin thickness 
of and the astronomical observations of the radius of 0.5Mq neutron star, is also suggested 1 ^ . In this paper, we adopt 
the normalized coupling constant such as SttkEs, where Sg denotes the saturation density, i.e.. Eg = 2.68 x 10^^ g/cm^. We 
remark that SttkEs becomes a dimensionless parameter. 

EiBI is characterized by two independent fields, i.e., the physical metric and the connection E^^. So, varying the action 
with respect to E^^ and one can obtain the held equations for A = 1; 

Qfii> — 

where q is determinant of and q^i, is an auxiliary metric associated with the connection as = 

{qaa,p + qa/ 3 ,a — qap,!j) /2. denotes the energy-momentum tensor, which is given by = {5Sm/ 5g^,/)/—g. 
Equation (O shows that the auxiliary metric q^^i, becomes equivalent to the physical metric if = 0. That is, EiBI with¬ 
out matter reduces to general relativity in vacuum 1^ . In addition to the above held equations, the energy-momentum tensor 
should satisfy the conservation law, i.e., = 0, where the covariant derivative is defined by g^^. As far as we know, 

unfortunately, there is no explicit proof that the conservation law of V = 0 is directly derived from the held equations (| 2 |l 
and dp. However, since the matter held is minimally coupled to the metric the conservation law might be obtained as in 
Ref. IIT^ . if the argument in ll3^ is applicable even for a bi-metric theory like EiBI. 

Now, we consider the neutron star models in EiBI. In general, magnetized neutron stars could deform due to non-spherically 
symmetric magnetic pressure. However, the magnetic energy in the neutron star is much smaller than the gravitational binding 
energy even for a magnetar, which is strongly magnetized neutron star. That is, the deformation due to the magnetic pressure is 
quite small and the shape of star is almost spherically symmetric. Thus, in this paper, we neglect the stellar deformation induced 
by the existence of magnetic held. Under such assumption, the equilibrium stellar model can be determined as a solution of 
Tolman-Oppenheimer-Volkov (TOV) equations in EiBI ll25l - l^ . The metric describing the stellar models is given by 

g^^dx^dx'' = -h f{r)dn^, (4) 

q^.^dx^^dx'^ = (5) 

where = dO"^ +sin^ 9d(p. In this paper, we consider the stellar models composed of perfect fluid, i.e., = {e+p)u^u'' + 

pg^'', where e, p, and are the energy density, pressure, and four velocity of matter given by 0, 0,0). Then, one 


( 2 ) 

(3) 
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FIG. 1; Neutron star models in EiBI with FPS EOS. The left panel corresponds to the stellar mass as a function of the central density normalized 
by the saturation density, while the right panel corresponds to the stellar mass as a function of the stellar radius. The solid line denotes the 
result in general relativity (k = 0) and the other lines denote the results in EiBI with various normalized coupling constant SttkEs- 


can show that abf = from Eq. where a and b are given by a = v/r+"8¥^ and b = \/l — SnKp, respectively 1^ . In 
addition to TOV equations, one needs to prepare the EOS for neutron star matter to construct the stellar models. We particularly 
adopt EPS and SLy4 EOSs ll^ . which are based on the Skyrme-type effective interaction (also see for the adopted 
EOSs). Eigure[T] shows the neutron star models constructed with FPS EOS, where the left panel corresponds to the stellar mass 
as a function of the central density normalized by Sg, while the right panel corresponds to the stellar mass as a function of the 
stellar radius. In this figure, the solid line denotes the results in general relativity and the other lines denote the results in EiBI 
with various values of SttkEs- From this figure, one can easily observe that the mass and radius of neutron stars depend strongly 
on the coupling constant in EiBI, even if EOS of neutron star matter is fixed. In fact, the stellar radii of IAMq neutron stars in 
EiBI become 9.3% smaller for SttkSs = —0.02, 7.6% larger for SnuSs = 0.02, and 16.5% larger for SttkEs = 0.05, compared 
with that in general relativity. 

On such a neutron star model, we consider an axisymmetric magnetic field generated by a four currency J^, adopting an ideal 
MHD approximation. The electromagnetic field is governed by the Maxwell equations with the physical metric 

= 0 , ( 6 ) 

(7) 

where is the Faraday tensor and the covariant derivative would be calculated with the physical metric Equation (|6ll 
automatically holds by introducing a vector potential, associated with as F^j^ = — A^ ,j. In order to determine 

the geometry of magnetic field, one also needs the equation of motions in addition to Eq. which is obtained by projecting 
= 0 on to the hypersurface normal to u^. With the ideal MHD approximation, the equation of motions becomes 

(s “F -\-11^11 p^jj = F^iiJ . (8) 

Now, assuming the appropriate gauge condition, Afj, can be described as = (0, Ar, 0, A^). In general, A^j, can be expanded, 
such as 


A^{r, 9) = ae{r) sin 9 dePi{cos 9), (9) 

where F^ (cos 9) is the Legendre polynomial of the order £. Furthermore, we especially focus on the dipole magnetic field, i.e., 
£ = 1, because the dipole fields could be dominant in the neutron stars. Then, in the same way as in Refs. IHdll^ . one can 
derive the equation to determine the vector potential oi; 

a" + - a\ + ^ e^oi = —47re^ji, (10) 

where the prime denotes partial derivative with respect to r, ji = co/(e + p), and cq is constant. We remark that constant C 
in Eq. (fTOl) is associated with the radial component of vector potential, i.e., Ar = . The procedure how to 

derive Eq. (fTOl l is detailed in Appendix lAl Consequently, since the magnetic field can be given by F°‘^/2, the 
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components of the magnetic field are expressed as 


Br =cos 9, (11) 

Bg = —sin0, (12) 

=-Caie"‘'/2 sin^ 0, (13) 

where e^i/a /3 denotes the totally antisymmetric tensor and etre<i> = yf—Q- From these expressions, one can see that the constant 
C, corresponds to the strength of toroidal magnetic field. Additionally, the tetrad components of magnetic files are given by 

i?p] = 2ai/“^cos0, (14) 

=-a'i/“^/^e“^/^sin6», (15) 

-S[0] = -Cai/“^^^e“'"/2sin6». (16) 


Since we consider that the exterior region of the star is in vacuum, as mentioned before, the spacetime outside the star 
becomes the same as that in general relativity, which can be described as the Schwarzschild metric. In such spacetime, the 
poloidal magnetic field (C = 0) is analytically given by 

(ex) 3/i^)r 

where /if, is the magnetic dipole moment observed at infinity ll43l] . Thus, at the stellar surface, the interior solution determined 
from Eq. ( fTOl i should be connected to the exterior solution [Eq. ( fTTl il in such a way that oi and a[ become continuous. In practice, 
from Eq. (O, one can show that the behavior of ui in the vicinity of the stellar center is expressed as oi = aor^ where 

ao is an arbitrary constant. So, the arbitrary constants ao and cq, which is a constant in the four currency ji, are determined so 
that oi and a'^ should be continuous at the stellar surface. 


In 1 - 


2M 


+ 


2M 2M^ 


(17) 


III. NUMERICAL RESULTS 

The magnetic field strength B is calculated by B , which can be expressed as 

B = f~^ [da^ cos^ 6 + sin^ 9 + CYa\fe~'^ sin^ 0] . (18) 

Thus, one can show that the magnetic field strength at the stellar center is Bq = 2aoao6o, where ao = a/1 + Sttkso and 
bo = a/1 — Sttkpo, while eg and po denote the central values of e and p. In the limit of k = 0, this expression reduces to that 
in general relativity It^] . The concrete structure of magnetic fields is discussed below, where we separately examine the pure 
poloidal magnetic fields (( = 0) in IIIII Al and the mixed magnetic fields (( Y 0) in ^IUBI 


A. Pure Poloidal Magnetic Eields (^ = 0) 

Eirst, one can show that the magnetic distribution is scaled by the magnetic field strength at the stellar surface of the poles 
(9 — 0), Bp, if the stellar model is fixed. That is, the distributions of B^^ /Bp for i = r, 9, and / are independent of Bp for each 
stellar model. In Eig. |2] we show the distributions of B[r]/Bp on the symmetry axis {9 — 0) in the left panel and B^g^Bp on 
the equatorial plane {9 = 7r/2) in the right panel for the stellar models with M — 1.4M0 contracted with EPS EOS, where the 
solid line corresponds to the result in general relativity and the other lines correspond to the results in EiBI with various values 
of SttkEs- Erom this figure, one can observe that the magnetic distributions in EiBI are qualitatively the same as that in general 
relativity. In fact, the deviation between the results in general relativity and in EiBI is not so much. In Eig. [2 we show the 
relative deviation of in EiBI from that in general relativity for the stellar models with M = IAMq constructed with 

EPS EOS. Erom this figure, we find that the deviation from general relativity is at most 10 % with the coupling constant in EiBI 
adopted in this paper. In particular, the magnetic distribution in the crust region depends weakly on the coupling constant in 
EiBI, which is less than 0.5%. That is, apart from the gravitational theory, one might be able to discuss the magnetic properties in 
the crust region of neutron stars. In addition, we show the magnetic configurations on the meridional plane for the stellar models 
with M = 1.4Mq for EPS EOS in Eig. |4] where the middle panel corresponds to that in GR (k = 0), while the left and right 
panels correspond to those in EiBI with Sttk^s = —0.02 and 0.05. The magnetic field strength is normalized by the magnetic 
dimple moment. As shown in Eig. |2] the magnetic configurations in EiBI are quite similar to that in GR. As with Eig. |2] we also 
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FIG. 2: For the stellar models with M = 1.4M0 for FPS EOS, the tetrad components of the pure poloidal magnetic fields are plotted as a 
function of r/R, where the left and right panels correspond to the radial component on the symmetry axis (6 = 0) and the 0-component on 
the equatorial plane (9 = 7r/2), respectively. The both components are normalized by Bp, which is the magnetic field strength at the stellar 
surface of the poles. The solid line corresponds to the result in general relativity, while the broken and dotted lines correspond to the results in 
EiBI with various values of SttkEs. The vertical lines denote the position of the stellar surface. 
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FIG. 3: Relative deviation of in EiBI from that in general relativity for the stellar models with M = IAMq constructed with FPS 

EOS, which is corresponding to the left panel in Fig.|2 The labels in the figure denote the values of coupling constant in EiBI. 


show the magnetic distributions for the stellar models with M = IAMq constructed with SLy4 EOS in Fig. |5] Comparing Fig. 
12] with Fig. |5] we find that the dependence of magnetic distribution on the coupling constant in EiBI is comparable to that on 
EOS for neutron star matter. 

Moreover, the magnetic field strength at the stellar center can be also scaled by Bp for each stellar model, such as 


Bo = PBp, (19) 

where /3 is a proportionality constant jidl] . In Fig. |6| we show the proportionality factor /3 as a function of the stellar mass with 
various values of the coupling constant in EiBI, where the left and right panels correspond to the results for the stellar models 
constructed with FPS and SFy4 EOSs, respectively. From this figure, one can see that the value of /3 is almost independent of 
the adopted EOS for neutron star matter, which is ^ 5. Additionally, for the stellar models whose masses are smaller than a 
critical value depending on the adopted EOS, /3 for the fixed stellar mass is almost proportional to the coupling constant in EiBI 
at least in the range adopted in this paper. This statement is clear from Fig. |7] where the proportionality factor in Eq. ( [19] ) for the 
fixed stellar mass is shown as a function of the coupling constant in EiBI. From this figure, such a critical stellar mass would be 
around 1.2Mq for FPS EOS and IAMq for SLy4 EOS. Furthermore, from this figure, we find that /3 for each coupling constant 
SttkEs depends weakly on the EOS for neutron star matter, if the mass of neutron star would be very low, for instance M ~ Mq. 


B. Mixed Magnetic Fields (C 7^ 0) 

As with the case of the pure poloidal magnetic fields shown in the previous subsection, the distribution of the mixed magnetic 
fields is also scaled by Bp, and the profiles of B[j] /Bp for i = r, 9, and (p are independent of the strength of Bp for each stellar 
model. For reference, first, we show the magnetic distributions in general relativity for the stellar models with M = \AMq 
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FIG. 4: Magnetic configurations on the meridional plane for the stellar models with M = 1.4M0 for FPS EOS in GR (middle panel) and in 
EiBI with SttkEs = —0.02 (left panel) and 0.05 (right panel). The magnetic field strength is normalized by the magnetic dipole moment 



EIG. 5: Similar to Eig. but for the stellar models constructed with SLy4 EOS. 


constructed with FPS EOS in Fig. |8] In this figure, the left, middle, and right panels correspond to B^r] /Bp on the symmetry axis 
{6 = 0), /Bp on the equatorial plane (9 = 7r/2), and /Bp on the equatorial plane, respectively. The solid line denotes the 
magnetic distribution for the pure poloidal field, while the other lines denote those for the mixed fields. As mentioned before, the 
toroidal magnetic component is characterized by the parameter (/, such as Eq. (fThl) . In Fig. 0 we show the magnetic distributions 
with the variable values of C normalized by 1/i?, because is a parameter with the dimension of inverse of length and then (/R 
becomes a dimensionless parameter. From this figure, one can observe that the distributions of B[r]/Bp and B^g-^/Bp are also 
changed due to the existence of the toroidal magnetic field, where the central field strengths of B[^]/Bp and B[g]/Bp decrease 
with the value of (. This result suggests the existence of the maximum of (/, where the central values of B[r]/Bp and B[g]/Bp 
become zero. Hereafter, such a maximum of ( denotes Cmax^ and Cmax^ = 3.30 for the case of the neutron star model in Fig. 

In practice, with ( more than Cmax, the direction of the magnetic field can be opposite inside the star ll4lll^ . Additionally, 
from Fig. [8] one can see that the position where | /i?p | becomes maximum is shifting outward with the value of (/. 

On the other hand, the magnetic distributions of B\^^-\/Bp, Bt^g^/Bp, and B^^-^/Bp for the stellar models in EiBI with various 
coupling constants are shown in Fig. |9l where the upper, middle, and lower panels correspond to the results with SttkEs = —0.02, 
0.02, and 0.05, respectively. Comparing this figure with Fig. one can see that the profiles of magnetic distributions in EiBI 
are basically similar to that in general relativity, where the distributions of B^^^/Bp and B^g^/Bp depend strongly on that of 
B[^-^/Bp. We also find that, as the coupling constant becomes smaller, the magnetic distributions are more sensitive to the value 
of (/R. For example, for (/R = 3, one can see that B[r]/Bp and \B[g-^/Bp\ with SttkEs = —0.02 in the vicinity of stellar center 
become smaller than that in general relativity. As a result, it is expected that the value of Cmax^ for the stellar model with smaller 
coupling constant in EiBI could be smaller. In fact, we find that Cmax^ = 3.07, 3.30, 3.41, and 3.51 for the stellar models with 
M = 1.4Mq constructed with FPS EOS with SttkEs = —0.02, 0, 0.02, and 0.05, respectively. 

Furthermore, in Fig. [Tol we show the values of CmaxT? for the stellar models with various stellar masses from M = Mq up 
to the maximum mass, where the left and right panels correspond to the results for the stellar models constructed with FPS and 
SFy4 EOSs, respectively, and the labels in the figure denote the values of the coupling constant in EiBI. From this figure, one can 
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FIG. 6: The proportionality factor P for the stellar models with various values of coupling constant in EiBI are shown as a function of the 
stellar mass, where the left and right panels correspond to the results for FPS and SLy4 EOSs, respectively. The labels in the figure denote the 
values of coupling constant in EiBI. 




FIG. 7: The proportionality factor /I as a function of the coupling constant in EiBI for the fixed stellar mass, where the left and right panels 
correspond to the results for FPS and SLy4 EOSs, respectively. 


see the value of Cmax^ for low-mass neutron star model is almost independent of not only the EOS but also the coupling constant 
in EiBI. On the other hand, the magnetic fields for the stellar models with canonical mass are more or less dependent on both 
the EOS and the coupling constant in EiBI. That is, an uncertainty due to the EOS for neutron star matter is degenerate into that 
due to the coupling constant in EiBI. Thus, only the measurement of the magnetic properties for the neutron star with canonical 
mass might be insufficient to observationally distinguish EiBI from general relativity. Anyway, the additional observations of 
the relativistic objects must become important to probe the gravitational theory in the strong field regime. 



FIG. 8: For the stellar models in general relativity with M = IAMq constructed with FPS EOS, the tetrad components of the mixed magnetic 
fields normalized by Bp are plotted as a function of r/R, where the left, middle, and right panels are B^^j/Bp on the symmetry axis (6 = 0), 
B[g]/Bp on the equatorial plane (6 = tt/2), and B[^]IBp on the equatorial plane, respectively. The different lines in the figure correspond to 
the magnetic field distributions with different values of C^R, and the labels in the figure denote the value of (R. 
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FIG. 9: Similar to Fig. but in EiBI with various coupling constants. Upper, middle, and lower panels correspond to the results in EiBI with 
StthEs = —0.02, 0.02, and 0.05, respectively. 



EIG. 10: Maximum value of allowed for the stellar models with various coupling constants in EiBI are shown as a function of the stellar 
mass, where the left and right panels correspond to the results for EPS and SLy4 EOSs, respectively. The labels in the figure denote the values 
of the coupling constant in EiBI. 


IV. CONCLUSION 

We consider the magnetic fields in the neutron stars in EiBI, where we especially focus on the dipole magnetic fields because 
such helds must be dominant in old neutron stars. To construct magnetic helds inside the neutron star, we derive the relativistic 
Grad-Shafranov equation in EiBI. Since the spacetime in vacuum in EiBI is equivalent to that in general relativity, i.e., the 
Schwarzschild spacetime, the magnetic held in EiBI outside the star is also equivalent to that in general relativity. In such a 
way that the interior magnetic helds should be connected to the exterior solution, the structure of magnetic helds is determined. 
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Then, we find that the magnetic geometry inside the neutron stars in EiBI is qualitatively similar to that in general relativity. The 
deviation of magnetic fields in EiBI from that in general relativity is not so much, which is almost comparable to the uncertainty 
due to the EOS for neutron star matter. Therefore, it might be difficult to distinguish EiBI from general relativity only by using 
the observations of the magnetic properties in neutron stars. However, the magnetic fields in the crust region for the neutron 
star with canonical mass depend weakly on the couplii^ constant in EiBI, while the crust properties such as the crust thickness 
depends strongly on the EOS for neutron star matter 0. That is, independently of the gravitational theory, one might be able 
to see the information about the EOS in crust region through the observations associated with the phenomena in crust region, 
such as the stellar oscillations. Anyway, there are many uncertainties in the magnetic properties even in general relativity, such 
as the magnetic geometry and the current distribution supporting the fields, although we consider only dipole magnetic fields in 
this paper. Comparing to such uncertainties, the imprint of EiBI gravity on the magnetic fields is weak, which suggests that the 
magnetic field could be a poor probe of gravitational theories. 
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Appendix A: Derivation of Eq. flOt 


According to Refs, llddlld^ . we briefly show in this appendix how to derive Eq. Coll, which is the equation to determine the 
distribution of magnetic field inside the star. One can obtain the following equations from Eq. Gli 


47r sin^ 9 


while from Eq. ® 


Anr = 


Ar 


99 


• Ar 




Ar,9r “ 1 “ A'p 


1 

-T 


,-A 


A, 


(/>,rr 


v' \' 
"2 ~ "2 


cos 9 
sin0 

^ _ y 

— -pA^,e0 


1 cos 9 

P‘2 • ^ A(p 0 , 

smO 


u' 

Ar,6r +A^,6J^ =0, 

41^,. J"+ 41^.0 J" = 0. 

Equation (IA6b with Eqs. (lAll i and (IA21 i can be written as 

9,0A(p^r A — 0 , 


(Al) 

(A2) 

(A3) 


(A4) 

(A5) 

(A6) 

(A7) 


where p = Arfi sin0. Thus, p should depend only on as rj = with a constant (. As a result, Ar is expressed 

as Ar = (e if A^ is expanded as Eq. ®. 

On the other hand, Eqs. (IA41 l and (IA5b become 


X,r — A^^rJ ^, 

where x J ^6 defined as 

, sV' , 

X,r = , 

J = - ' . 2 /I i J<t>~ ■ 

f{e+p)siir9 V 47r 
Owing to the relation x,r6 = X,0r with Eqs. dASb and (IA9b . one can obtain 

A(j)^rJ,0 — A^^sJ^r = 0 . 


(A8) 

(A9) 


(AlO) 

(All) 


(A12) 
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Therefore, J also depends only on A^, i.e., J = —cq — ciA^, where cq and ci are constants. That is, 

J^=^e ’'A^-{co + CiA^)f{e+p)sin^0. (A13) 

47r 

At last, substituting Eqs. (|9]) and (IA131 l into Eq. (IA31 l. one can obtain the equation (fTOl i describing the function of oi We remark 
that the term of ci in Eq. (IA13b is neglected, because we focus on the dipole (i = 1) magnetic distribution in this paper and the 
term of ci can contribute as the multipole higher than f = 3 lldHld^ . 
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